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Nucleon-nucleus optical potential in the particle-hole approach 
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Feshbach's projection formalism in the particle-hole model space leads to a microscopic description 
of scattering in terms of the many-body self-energy. To investigate the feasibility of this approach, 
an optical potential for 16 O is constructed starting from two previous calculations of the self-energy 
for this nucleus. The results reproduce the background phase shifts for positive parity waves and 
the resonances beyond the mean field. The latter can be computed microscopically for energies of 
astrophysical interest using Green's function theory. 

PACS numbers: 24.10.cn, 25.40.Cm, 25.40.Lw, 21.10.Jx, 21.60.Jz. 



I. INTRODUCTION 

The Feshbach's theory of the optical potential 1] pro- 
vides a tool to describe the scattering of nuclcons from 
nuclei. In its original form, the full Hilbert space is parti- 
tioned in a subspace that contains only one particle added 
to a particular state of the core nucleus. This includes 
the elastic scattering states and the nuclear orbits un- 
occupied by the target's nucleons. The resulting optical 
potential can be thought as an effective interaction that 
accounts for the effects of the degrees of freedom of the 
excluded space — such as the overall antisymmetrization 
of the wave function, excitations of the target or breakup 
channels. The complexity of these effects makes the ab 
initio calculation of the optical potential a very difficult 
task. The choice of working in the space of one parti- 
cle plus a core is shared by several theories of nucleon- 
nucleus scattering. Examples are the cluster model Q, 
folding potential jj] and shell model embedded in the con- 
tinuum 4]. All these techniques have been successfully 
applied to scattering processes. 

A conceptually different approach consists in extend- 
ing the scattering space to include both a particle on top 
of the nuclear core and the possibility of propagating a 
hole excitation. In the following, we will refer to this as 
the particle-hole (ph) Hilbert space. Mahaux and Ca- 
puzzi 5] and Jennings and Escher |(J have shown that 
applying the Feshbach's formalism to this space leads to 
an optical potential that is the usual many-body self- 
energy defined in Green's function theory 0, II]. The 
properties of using the self-energy as an optical potential 
have been discussed by Mahaux and Sartor in Ref. 0. 
The solutions of the scattering equation above (below) 
the Fermi energy are the overlap wave functions between 
the core and the eigenfunctions of the systems with A + 1 
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(A — 1) nucleons. For particle states these are the same 
wave functions as obtained in the original Feshbach ap- 
proach The fact that states with A — 1 particle are 
included in the formalism does not lead to complications 
in the study of scattering events. Rather, more physical 
information can be extracted using the self-energy since 
it also describes the so called 'Pauli forbidden' orbitals, 
occupied by the nuclcons of the target. Correspondingly, 
the optical potential in the ph space is more easily compa- 
rable to phenomenological models (based on, e.g., Wood- 
Saxon wells) 0, |i| since these also describe hole states. 
It also has better analytical properties than its counter- 
part in the particle-only space ja, |(| . Finally, the theory 
of Green's functions provides a natural way to include 
the effects of the excitations of the core in terms of an 
expansion in Feynman diagrams 0, [Toj . 

Many-body Green's functions have been applied in the 
past to study nuclear correlations, with emphasis on the 
hole part of the one-body spectral function. Recent de- 
velopments offer the opportunity to obtain sophisticate 
descriptions of the couplings between particles and collec- 
tive states using a Faddeev expansion llulTi . A similar 
formalism was already considered in Ref. [13j for the opti- 
cal potential where, however, the dressing of propagators 
was disregarded and no application was attempted. The 
calculations of Refs. [Til Il4| account for collective mo- 
tion near the Fermi level, including the energy regime of 
interest to nuclear astrophysics. With a variety of exotic 
isotopes that are involved in stellar processes becoming 
experimentally accessible in modern radioactive beam fa- 
cilities, it is important to investigate whether the Green's 
function approach can be applied to study low-energy nu- 
cleon scattering and capture processes. 

We have recently considered the self-energy resulting 
from a recent application of the self-consistent Green's 
function (SCGF) method to 16 O [3 and explored its 
predictions for proton-nucleus scattering. These results 
were obtained in a restricted model space, which is not 
fully appropriate to describe scattering events. However, 





FIG. 1: Feynman diagrams representation of the self energy. 
The first diagram on the r.h.s. represents the Hartree-Fock 
like contribution to the mean field. The remaining ones de- 
scribe core polarization effects in the particle (2pfh) and hole 
(2hf p) part of the spectrum. 



preliminary calculations with only this input gave en- 
couraging results 0| . In this work the self-energy of 
Ref. |l4| is augmented by including the components from 
outside this model space, as they have been computed in 
Ref. [16j ■ We then report on the final conclusions of these 
exploratory studies. 

The model and the details of the calculations are given 
in Sec.|n] The results for elastic scattering and the bound 
states of 17 F are reported in Sec. IIIII and are preceded 
by a discussion of the different contributions to the self- 
energy, in Sec. lIII K\ A discussion and our conclusion are 
given in Sees. II VI andlVl 



II. MODEL 

In the SCGF approach it is useful to split the many- 
body self-energy into three contributions [Tol as 
shown diagrammatically in Fig. ^ There, the dou- 
ble lines represent the exact one-body Green's function, 
which contains complete information on the particle and 
hole spectral distributions. The first diagram on the 
r.h.s. is the direct extension of the Hartree-Fock poten- 
tial to include the effects of the fragmentation of strength 
and represents the nuclear mean-field (MF) in the pres- 
ence of correlations. The remaining contributions split 
naturally in diagrams containing at least two-particle- 
one-hole (2plh), describing the system of A+l parti- 
cles, or two- hole-one-particle (2hlp), corresponding to 
A-l particles. The irreducible propagators i?( 2 ? lft ) and 
R (2hi P ) account f or tn e core polarization contributions 
to the opticalpotential in the particle and hole spaces, 
respective ly El . The separation of Fig. ^ is exact. In 
Refs. jTl 137 R {2plh) and Ri 2hl P) were computed em- 
ploying a Faddeev expansion that permits the direct cou- 
pling of the single-particle motion to collective excita- 
tions of the core. These were evaluated in the dressed 
random phase approximation (DRPA) Jl£|. The exam- 
ple of a diagram that contributes to R\ 2 P xh > is given in 
Fig. |21 Since this expansion is based on the fully frag- 
mented single particle propagator — which is generated 
from the self-energy itself — a self-consistent solution is 
required. 

The nuclear self-energy computed in Ref. |14| was ob- 
tained within a model space V corresponding to the har- 




FIG. 2: Example of a diagrammatic contribution included in 
the Faddeev expansion for R ( - 2 P lh 1 (see Fig.[TJl. A quasiparti- 
cle is coupled to the response function 'n} ph > that describes the 
target nucleus. It can also participate in a pairing processes, 
which is accounted by the two-body propagator g u ^ pp \ 



monic oscillator wave functions for all orbitals up to the 
pf shell plus the g 9 / 2 - A parameter 6=1.76 fm was em- 
ployed. This space appears to be large enough to describe 
the influence of the low energy (long-range) excitations 
on nuclear fragmentation |12| . However, it requires a 
proper extension for applications to single particle scat- 
tering, as it will be discussed below. The effect of cor- 
relations outside this model space were accounted for by 
employing a G-matrix as an effective interaction, which 
was derived from the Bonn-C potential |19|| according to 
Ref. [2(J. The computation of the G-matrix for posi- 
tive energies is an outstanding problem which was not 
attempted there. Therefore, we employed a fixed start- 
ing energy of -5MeV in the present work, as the closest 
reliable choice to the continuum. 

At low energies the optical potential is well approxi- 
mated by a real interaction and i?( 2 P lh ) and R( 2hl P} can 
be expressed as discrete sums of poles. Thus, for each 
given partial wave, Ij, the contributions depicted in Fig.^ 
can be expressed as 



^MF,Fadd 
J l3 



(k,k')= e 



(la) 



^{2plh),Fadd^ ^ 



E 



n a -na £V 



E 



(ib) 



^(2hlp),Fadd^ ^ 



E 



E 



(K y 



- u-e u -IT] 



(lc) 



3 



where 4> a {r) are the harmonic oscillator radial func- 
tions referring to sing le particle quantum numbers 
a = {n a , l a ,j a , ma] |38j . the first sum runs over all the 
orbits belonging to the model space and l a j a — Ipjp = Ij 
since 16 has a + isoscalar ground state. 

The superscript "Fadcf indicates that Eqs. repre- 
sent the results of Ref. This is the most sophisti- 
cated calculation available to date for the self-energy at 
low-energies that account for the coupling between single 
nucleons and collective excitations. However, the expan- 
sion over a few harmonic oscillator states is not optimal 
for describing the details of the nuclear surface. Anal- 
ogously, it misses part of the large momentum compo- 
nents in the optical potential. This is particularly critical 
for the MF component, which describes the background 
of the phase shifts. On the contrary, the same nucleus 
was studied in Refs. |16| employing a spherical box ba- 
sis that includes all the relevant momentum components. 
An effective G-matrix, derived for nuclear matter and 
the Bonn-B potential [l!j , accounted for the binding due 
to short-range and tensor correlations. The self-energy, 
computed only to second order in the perturbation series, 
neglected most of the collective effects. This approach 
was applied to obtain the quasihole wave functions asso- 
ciated to the p states occupied in 16 O, with sufficiently 
accurate results to describe the shapes of the (e, e'p) cross 
sections to those states pi) . 

In this work, we chose to employ a mixed represen- 
tation of the self-energy in which the MF components 
missing in the space V were extracted from Refs. Il6l. 
while the contributions beyond MF computed in Ref. |l4f 
[Eqs. (|Tb|) and (fTc| j_ were retained. To do this the MF 
self-energy of Ref. [lfj was split in two parts, 
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where E, 

acts on the excluded space. Two approximations were 
considered depending on which MF component to employ 
inside V. In the first case (I), ^ IF ' Box was added to 
Eq. (|lafl . In doing this, we note that the G-matrix used 
to compute Yi MF ' Fadd accounts for the extra binding due 
to the degrees of freedom of the excluded space. Since 
these are reinserted explicitly by Y^^ F ' Box , one should 
also rescale Y, MF ' Fadd appropriately by a constant, N 1 . 
The second choice (II), consisted in employing both parts 
of Eq. (0). Also in this case we kept the possibility of 
tuning the depth of the potential. The complete MF 
contributions employed in this work are 
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where the constants Ay and ATy depend of the specific 
channel and will be discussed below. The full self-energy 
employed in the calculations is [see Fig. ^ 

= Ey*™(M) 
+T l ^ plh) ^ Fad \k,k l -u J ) + T l ^ hlp ^ Fadd {k,k'-u J ) . (4) 



The Dyson equation can be expressed in a Schrodinger- 
like form, where the self-energy takes the place of a non- 
local and energy dependent optical potential [h — c = 1 
and /i is the reduced mass] 
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where Vq ou1 (k, k') in the Coulomb interaction cor- 
responding to a uniformly charged sphere of radius 
R c = 3.1 fm. This was added to account for the elec- 
tromagnetic interaction missing in the calculations of 
Refs. Ill E3. Due to the non local character of E* 
Eq. JSJ is conveniently solved in momentum space. In 
doing this, the long distance part of the Coulomb po- 
tential was solved using the Kwon-Tabakin-Lande |23| 
procedure for bound states and the Vincent-Phatak p3l | 
one for scattering. 



Above the Fermi level the eigenvalues of Eq. (|SJ) are 
related to the spectrum of 17 F by E" m = E n F — E g s °. 
Thus, E cm > describes the scattering of protons from 



16 O while the bound solutions are the overlaps of the 
ground state of 16 O with the corresponding bound states 
17 F. Analogously, below the Fermi level E™ n 



E^ N and the eigenstates represent the overlaps with 15 N. 
The Dyson equation implies that the bound solutions of 
Eq. ||SJ have to be normalized to their spectroscopic fac- 
tor according to 
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where ip n {k) is the solution itself normalized to unity 
and E" is the corresponding eigenvalue. The asympotic 
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E%* {MeV) 


Si /o 


-3.57 


0.69 


1.05 


-0.1 


^3/2 


1.87 


0.72 


1.08 


4.4 


P3/2 


-16.61 


1.06 


1.07 


-18.5 



0.95 (£^ /2 =-15.1) 



a % = l ' e . 

b E l J m = Ecm, except when specified. 

TABLE I: Corrections applied to the depth of the MF poten- 
tial E MF ' J [Eq. ga)] and quasi particle energies obtained in 
the calculations of Fig. 



normalization for the unbound solutions is related in the 
usual way to the flux of incoming particles. 



III. RESULTS 

Eqs. (|3J and (0J include the relevant physics from both 
the calculations of Refs. 0] and [l^. This self-energy 
represents a model for the optical potential that acts on 
the full ph Hilbert space and can give sensible predictions 
near the Fermi level. However, the two-body realistic in- 
teractions alone, as used in these works, cannot repro- 
duce the experimental bindin g en er g ies and spin-orbit 
splitting for nuclei with A > 3 |24tl25| . To obtain these, 
relativistic effects or three-body forces are required [26| . 
In this work £* was constrained to reproduce the exper- 
imental spectrum in two ways. First, the constants 

and N[j that affect the depth of the optical potential 
were chosen to reproduce the corresponding quasiparti- 
cle energies. These are the Si/ 2 and d 5 / 2 bound states of 
17 F, its dz/2 resonance and the px/z and P3/2 hole states 
of 15 N. Second, complex resonances that do not have a 
mean field character are generated by the dynamic part 
of the self-energy. At low energy, most of these couple 
to only one pole e l± in Eqs. (|lb|) and (|Tc|l . Therefore, 
we have fitted those poles that could be identified with 
specific resonances of the A+l system ( 17 F) by impos- 
ing that Eq. JSJ yields the corresponding experimental 
energies. We note that a similar approach was already 
employed in Ref. 0] . This is necessary for the particular 
case of 16 O due to the strong coupling between the single 
particle spectrum and collective motions, which suggest 
the need for an improved description of the low-energy 
structure of this nucleus and more attractive effec- 
tive interactions [2?! ■ Although, satisfactory results can 
already be obtained in similar calculations for heavier 
nuclei [HIUIIl 

The influence of this fitting procedure on the results 
is discussed in the following. After calibrating Eqs. (JTJ 
and © to the spectra of 17 F and 15 N, the results for 
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FIG. 3: Phase shifts obtained from the self consistent self- 
energy of Ref. |3 [Eq. Q before (dotted lines) and after 



(dashed lines) correcting the depth of E 



to reproduce 



the quasiparticle (and quasihole) energies. The full lines are 
obtained when also the momentum components outside the 
model space V are included, Eq. 13all . The S1/2, ^3/2 an d 
P3/2 partial waves are shown. For P3/2, the dot-dashed line 
was obtained by fitting AjJ to reproduce the background 
phase shifts rater than the quasihole energy. The values of 
the corrections A; 7 are reported in Tab. CT] The experimental 
results are from Refs. |3ll | (crosses) and |32l | (circles). 



the scattering phase shifts and the bound single particle 
wave functions are a prediction of the model. 



A. Parameter dependence 

To discuss the influence of the different contributions 
to Eq. ((3J), the phase shifts for proton scattering have 
been computed employing different truncations of the 
mean field self-energy Y 1 MF - 1 . The results are shown in 
Fig. 13 for three partial waves. The dotted lines were 
obtained by retaining only the original contribution to 
the self-energy of Ref. [l4|. Thus, neglecting Y, 1 vrF ' Box in 
Eq. l(3"al) and setting Nh = 1 for all cases. The results 
obtained by constraining these constants to generate the 
proper quasiparticle energies is given by the dashed lines. 
The full line shows the full results form Eq. (|3af> . ob- 
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tained by including also the T, 1 



MF,Bc 



term and refitting 



the N 1 . The values for the quasiparticle energies and the 
constants used are given in Table [I] 

The background contribution to the phase shifts of the 
Sx/2 partial wave is described correctly by Y, MF ' Fadd but 
not the energy of the bound state. Vice versa, it is possi- 
ble to constrain the depth of the potential to reproduce 
the latter but the agreement with the experimental phase 
shifts is lost. However, both quantities are reproduced if 
YMF,Box . g mcm( j ec } j n ^ mg case t ne correction required 
in the depth of the potential, Af| = 1.05, is less sig- 
nificant than when only E MF > Fadd [ s included. A similar 
trend is seen for the <i 3 / 2 channel. Reproducing the en- 
ergy of the single particle resonance with Y l MF > Fadd alone 
requires a sizable change in its depth, while the observed 
phase shifts are obtained only after including the com- 
ponents outside the space V. We observe that the ex- 
pansion of Eq. includes only one harmonic oscillator 
function for the d 3 / 2 wave and two for S\i 2 - With such a 



restricted space, it is remarkable that the resulting back- 
ground phase shifts are still obtained somewhat close to 
the experiment. 

A different behavior is found for the I = 1 partial 
waves. The results for p 3 / 2 are shown in Fig. [3] (the P\i 2 
case is analogous). In this cases Y 1 MF ' Fadd produces a 
spurious resonance at ~1 MeV that is not seen experi- 
mentally. Fitting the potential's depth to constrain the 
quasihole energies of 15 N generates a more attractive well, 
thus worsening the situation. The phase shifts improve 
upon introducing Y^ IF ' Box (full line) but still show a rise 
of the background with the cm energy, while the experi- 
mental results are practically constant. A proper choice 
of JVi (and AfJ ) allows to reproduce the behavior of 

f»3/2 v Pl/2' r 

the phase shifts at the lower energies but results in un- 
derbinding the corresponding orbitals in 16 O (see Tab.Pl. 

The curves of Fig. |3| have been computed without any 
shift of the e l+ poles in in Eq. Ijlbfl . This gives an idea 
of the quality the energy spectra obtained adopting the 
interaction of Ref. [2(j • No solutions were obtained that 
could be interpreted as the d 3 / 2 resonances above 5 MeV. 
The S1/2 resonance at ~6 MeV was obtained as a cou- 
pling of a proton to the first excited state of 16 O. Anal- 
ogously, the two lowest resonances in both p 1 / 2 and p 3 / 2 
can be interpreted as quasiparticle interacting with the 
first isoscalar 3— and 1— levels of 16 1331 . 
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FIG. 4: Phase shifts for positive parity waves obtained from 
the the self-energies I (dashed lines) and II (full lines). The 
experimental results are from Refs. |3lt ) (crosses) and [32 ] (cir- 
cles). 



tive parity waves and it is more accurate for the P\/ 2 
case, for which the collective resonances are sharper. In 
general, the non MF resonances were predicted narrower 
than the experiment. This is probably related to the lack 
of momentum components outside the model space V in 
Eqs. (jlb|) and ljTc|) . which were not corrected as for the 
MF part of the self-energy. 

The values of Nf- and Nf? show that much smaller 
modifications are needed to force Yj*' 1 to reproduce 
the quasiparticle energies. This is consistent with the 
more sophisticate treatment of long-range correlations 
achieved in Ref. 0] . It is worth noting that the inability 
to describe both the bound energies and the scattering 
for the negative parity waves is consistent with the lack of 
three body forces in the present model, which are needed 
to reproduce their spin-orbit splitting. On the other hand 
the results with the MF potential II reproduce reasonably 
well both the quasihole energies and the phase shifts. 



B. Phase shifts for proton scattering 

Figures 0] and \5\ compare the phase shifts obtained 
from both the Self-energies I and II after constraining 
the quasiparticle energies and resonances to their experi- 
mental values. Table UTI shows the values of the constants 
Nfj and Nf? used to obtain these results. 

For the positive parity waves the background phase 
shifts are described equally well by both optical poten- 
tials. The potential T,*' 11 can also describe the nega- 



C. Bound overlap wave functions 

The overlap wave functions associated to the bound 
states of 17 F are shown in Fig. [S] for the two choices 
of Eq. P|l. The asymptotic behavior in presence of a 
Coulomb field is given by 



W_^,/+i/ 2 (r) 



(7) 



6 



-a 

'-H 



10 

9 

8 

7 

6 

5 

4 

3 
Z 

to io 




I p 3/2 partial wave 




3 4 5 

E [MeV] 

cm 



FIG. 5: Phase shifts for negative parity waves obtained from 
the the self-energies I (dashed lines) and II (full lines). The 
experimental results are from Refs. {31j (crosses) and |32| (cir- 
cles). 
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(MeV) 
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1.24 


4.4 
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1.05 


1.24 


-0.1 


^5/2 


1.16 


1.29 


-0.6 


Pl/2 
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1.15 


-12.1 


P3/2 


1.07 


1.25 
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a E l J m = -BcTrfi except when specified. 

TABLE II: Corrections applied to the depth of the MF po- 
tentials I and II [Eq. ©] for the calculation of Figs. 0]and|K| 



where W- v ^ + i/2 is a Whittaker function, r\ the Som- 
merfield parameter and Cy the asymptotic normaliza- 
tion constant (ANC). The spectroscopic factors, ANCs 
and root-mean-square radii obtained are given Tab. IIIII 
The Self-energy X*' 7 predicts larger radii and ANCs 
than X*' 77 , which pulls these orbitals more strongly in- 
side the nucleus. At the same time both choices yield the 
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FIG. 6: Radial part of the overlap wave functions between 
16 and the bound d 5 / 2 and s 1 / 2 states of 17 F. 
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(4 /2 ) 1/2 


£*' 7 


0.931 


-82.5 


5.86 


0.913 


1.07 


4.01 




0.921 


-73.9 


5.55 


0.909 


0.81 


3.70 



TABLE III: Spectroscopic factors, ANCs (in fm" 1/2 ) and 
root-mean-square radii (in fm) for the bound d 5 / 2 and Si/ 2 
orbitals of 17 F. 



same spectroscopic factors, implying equal occupancies. 
The depletion of these orbits is driven by the coupling to 
long-range collective excitations contained in ^ 2plh ^ Fadd 
and YS 2hlp '' Fadd . An additional quenching is expected 
from short-range and tensor correlations and was not 
accounted for in this work. This has been seen to be 
of about 10% for bound orbitals of several closed shell 
nuclei ^(j- However, what the strength of this reduc- 
tion should be for loosely bound nucleons, which can be 
largely localized at radii outside the nuclear surface, has 
not yet been investigated. 



IV. DISCUSSION 

In order to include the principals physics ingredients, 
an improved self-energy was constructed from those ob- 
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tained by two different calculations of the nucleus of 16 O. 
The effects of coupling a nucleon to collective modes were 
studied in Ref. 14] but in a restricted model space. The 
momentum components of the mean field potential out- 
side this space were instead extracted from Ref. [I^. 
Since these calculations are based on realistic two-body 
inter-nucleon reactions, the energy spectra cannot be ac- 
curately reproduced. Therefore, the model has to be con- 
strained phenomenologically to reproduce the experimen- 
tal spectra of the nuclei with A±l nuclcons. Resolving 
this situation may require the use of multi-nucleon forces 
and more appropriate effective interactions. 

The present results show that both the inclusion of 
all momentum components of the particle-hole Hilbcrt 
space and a proper treatment of long-range correlations 
are important to correctly reproduce the mean field op- 
tical potential. The coupling of single particle strength 
to long-range excitations is also responsible for the cre- 
ation of non mean field resonances. After constraining 
the prediction for the single particle energies, the phase 
shifts for the scattering of protons from 16 O were ob- 
tained in fair agreement with the experimental data, ex- 
cept for the background behavior of the of the p waves. 
The difficulties for these waves are accompanied by the 
issue of explaining the hole spectroscopic factors with 
the same parity extracted from (e, e'p) experiments |34) . 
In Ref. [14j the latter were linked to the particularly 
complicate structure of the low-energy spectrum of 16 O 
and further studies along this line have been initiated in 
Ref. |12(. It is plausible that the required improvements 
will resolve both the problems of spectroscopic factors 
and scattering phase shifts. Similar issues are expected to 
be beyond the requirements for reproducing most heavier 
closed shell nuclei (Mill. 



V. CONCLUSIONS 



This work investigates the possibility of describing 
nucleon-nucleus scattering employing the many-body 
self-energy as an optical potential. This corresponds to 
applying the Feshbach projection formalism to an Hilbcrt 
space containing both particle and hole states. 

The present results are a first attempt at computing 
scattering processes using the many-body Green's func- 
tions and required to introduce a certain amount of phe- 
nomenological corrections. However, it is shown that pre- 
dictions for the scattering of nucleons can be obtained 
working also in the particle-hole space. The present re- 
sults also give insight into the developments that will be 
needed to pursue reliable microscopic calculations of the 
optical potential. We feel that the overall quality of the 
results can be com para ble to other methods applicable 
at low energies |3ti l37[ when the role of the missing in- 
gredients, such as short-range correlations, is included. 
Thus SCGF could be considered as a valid candidate for 
the study of selected reactions at astrophysical energies. 
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